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Abstract: 

In this work, we employ the Bethe-Salpeter (B-S) equation to investi- 
gate the spectra of free diquarks and their B-S wave functions. We find 
that the B-S approach can be consistently applied to study the diqaurks 
with two heavy quarks or one heavy and one light quarks, but for two 
light-quark systems, the results are not reliable. There are a few free 
parameters in the whole scenario which can only be fixed phenomeno- 
logically. Thus, to determine them, one has to study baryons which are 
composed of quarks and diquarks. 



1 Introduction 

The subject of diquarks has attracted attentions of theorists of high energy physics for 
decades. The reason is obvious. Since baryons are composed of three valence quarks, the 
three-body system is much more complicated than mesons which are two-body systems 
of quark and antiquark. If the diquark picture is applied, namely two of the constituent 
quarks constitute a color-anti-triplet sub-system, the three-body problem can be reduced 
to a much simpler two-body system. However, one would ask whether QCD, which is 
responsible for interaction between quarks or antiquarks, favors such diquark structure of 
two constituent quarks. 

Recently, the topic on diquarks revives because it may bring up some direct phe- 
nomenological consequences. At relativistic heavy ion collisions (RHIC) with high temper- 



1 



ature and density, diquark production and even diquark condensation become hot topics 
[1]. Some authors suggest that pentaquarks are of a diquark-diquark-anti-quark structure 
[2, 3, 4, 5], and Zou et al. proposed that there is a pentaquark component in nucleons [6]. 
All the subjects concern a dynamics which results in the substantial diquark structure. 

In 1964 Gell-Mann first proposed feasibility of the diquark structure [7], then Ida, 
Kabayashi[8], Tassie et al.[9] applied the concept to study baryons. Later many au- 
thors carefully analyzed the quantum numbers based on the group theory [10, 11, 12]. 
Anselmino[13] et al. indicated that two heavy quarks or one heavy and one light quarks 
may constitute stable scalar or axial vector diquark of color anti-triplet [3, 4]. However, 
for the two light quarks, it is not very clear if a substantial diquark can exist. 

In this work we employ the Bethe-Salpeter (B-S) equation to study the spectra of free 
diqaurks which are composed of two heavy quarks or one-heavy-one-light quarks, we also 
try to extend this approach to study the diquarks of two light quarks. Our results show 
that just as for the light pseudoscalar pions, this framework does not work well, unless 
one considers extra contributions to quark masses. Dai et al. [14] used the method to 
study the pion structure and recently, Wang et al. [15, 16] employed the same approach 
to study the diquarks which contain only two light quarks. We will come back to this 
issue in the last section. While numerically solving the B-S equation, we have employed 
the method developed by Chang et al. [17], which is proved to be powerful and efficient 
in the numerical computations. 

This work is organized as follows. After this short introduction, we derive our formula- 
tions. In Sec. Ill, we present the numerical results, and then in the last section, we briefly 
discuss our results. 

2 Formulation 

In this work, we use the B-S equations to study the diquark structure. First we construct 
the Green's function of the diquark in the 4x4 matrix form and derive the corresponding 
B-S equation. 

2.1 Structures of diquarks 

One can define the B-S wave function of a diquark as [18, 5] 



where i, j, k are color indices which will be omitted later without misunderstanding. A 
Fourier transformation brings it into the momentum space as 



xUxi,X2) = % < 0|T(Vi(xi)V,^(x2))|L' >= e 



-iP-X^k 



(1) 




(2) 



where 



Ip'^ix) = CllF , C = i'jQj2,X = aiXi + 02^2, X = Xi - X2, 
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mi m2 , . 

ai = ■ ,02 = : ,P = pi+p2,q = a2Pi - Q;iP2- (3) 

nil + "1-2 "T-l + "1-2 

Thus the B-S wavefunction in the momentum space reads 

i^i - mi)xp{qW2 +m2)=J d^kViP, q; k)xp{k), (4) 

where V{P, q; k) is the kernel. Although the B-S equation has the same form as that 
for mesons, the kernels for the two cases are different. For the diquarks there are no 
annihilation diagrams which exist for meson case. In fact, the definition of the diquark 
B-S wavefuntion automatically eliminates annihilation diagrams"^. Greiner [19] indicates 
that the B-S equations for quark-quark system (diquark) and for quark-antiquark system 
(meson) have a formal symmetry. 

According to the method given in Ref. [20], one can obtain an equation group which 
contains four independent equations. 



(M - uip - uJ2p)^V{q^pJ = Atp{q'^JMQk)^tpi<lpJ , (5) 

(M + uip + u;2p)'Pp~{qpJ = -Kpi<lpJM<lpJ^2pi<lpJ > (6) 

and 

<^r(9^j=o, (7) 

<^;+(g^J = 0. (8) 
The B-S wavefunctions are normalized as following: 

The kernel for the B-S equation for diquark has been widely discussed by many 
authors[21] 

2 a 

I{r) = Vs{r) + Vb + 7'' ® iM^) = /3Ar + Vb - 7'' ® 7mo — • (10) 

To avoid the infrared divergence, one can introduce a convergence factor e""''. Thus the 
potentials are respectively^ 

K(r) = ^(l-e— •), (11) 

and 

9 n/ 

K(r) = -3^e— . (12) 

""^The four-point Green's function for diquark is defined as < Q\T{ip{xi)tl)'^(x2)^'^{x3)ip{x4)\Q > whreas 
for meson it is < 0|r(-(/;(2;i)V'(a;2)^(a;3)i/'(a;4)|0 >, since i>{x) and ^p''{y) cannot contract, the annihilation 
diagrams do not exist in the diquar case. 

^In fact, for the Coulomb term, there is no infrared divergence, but the hnear confinement, which is in 
the form of in the momentum space, diverges. Introduction of a converging factor is necessary and 
is benign for the Coulomb term. 
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In the momentume space the potential reads 



i^(q) = IC5/K(q)+7^®7MK(q) (13) 

where 

V.(q) = -(^^ + W»(q) + ^(5^. (14) 

and 

The running couphng constant as{q) is 

as(q) = ^2 (16) 

where a is a constant which freezes the running couphng constant at low momentum. 

Because diquark exists in the color-3 state, < 3|A'^A'*|3 >= ^ < OjA^A^jO >, where 
|0 > is the color singlet. Thus the coefficient of the Coulomb term, which originates from 
the onc-gluon exchange, is — ^ for diquarks, whereas it is — ^ for mesons. The linear 
potential comes from the non-perturbative effects of QCD, it might have different form 
for diquark and meson, and so far we cannot determine its exact value from the QCD 
theory. Thus we introduce a phenomenological parameter 13. Indeed, some authors[15] 
argued that for diquarks, the linear confinement might not exist at all, if so, /3 = 0. On 
other aspect, if the dynamics for the bound state of diquark is similar to that for mesons, 
one should expect (3 = 0.5, in analog to the coefficient ratio in front of the Coulomb term. 
We set the parameter to be free within a range of < /3 < 1.0. 

2.2 Scalar diquark system 

We carry out all the computations in the center-of-mass frame of the diquark, thus P = 
(M, 0) where P is the total four-momentum of the diquark and M is its mass. 

For 0+ diquark, the general form of its wavefunction is 

</f(q) = 7o?'i(q) + ?'2(q) + ^±b3(q) + 7o^±64(q)- (17) 

Here we define q± as (0, q) which is perpendicular to P. Substituting the wavefunction 
into the equations (7) and (8), we can obtain the constraint conditions 

[ui + u;2)(q — mim2 — UI1U2) 

and 

hM) = q!("^^+"^^)^3(q) ^^g) 

q^ — mim2 — iO\U2 

With these constraints, the B-S wavefunction of a 0+ diquark has the following form 

q^(a;i -a;2)^7o64(q) , c^{mi + m2)h^{<i) , , , / ^ , i u t \ 

^(q) = 7 — ^+^"^ ^^^+^±^3(q)+7o^±^4(q)- (20) 

(mi — m2j(q — mim2 — 0J1U2) q — m\m2 — UJ1LO2 
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Substituting them into eqs.(5) and (6), we obtain an equation group for the component 
functions: 

2q2[(q2 - mim2 - a;ia;2)64(q) - (m2Wi + mia;2)63(q)] 

UVl — LJi — Ul2) S 

— m\m2 — 0J\0J2 

(a;iu;2)(k2 - mim2 - ujiki^2k) 

+(ys - 2Vy){miUJ2 + m2a;i)(k^ - mim2 - uJikUJ2k)fl ■ k] 

+b4{k)[{Vs - 2Va){uJik - uj2k fq^'k^ + K(q^ + mi7n2 + u;ia;2)(k^ - mim2 - wifca;2fc)q • k]| , 

(21) 

2q2[(q2 - mim2 - wiW2)fe4(q) + ("i2'^i + mia;2)63(q)] 
{M + wi + a;2 j 2 



■ / 7 ^,1.2 '^^'^ ^ {himVs + 4K)(mi + m2){coi + u;2Wk' 

J (L0lU;2)(k^ - 17111712 - UJikLi>2k) I ^ ' 



q^ — 17111712 — CO1UJ2 



= 1 1 ^I^2 V {-HmVs + 4F.)(mi + m2){u, + U2)c^^k^ 

J (a;ic<;2)(k^ - mim2 - t^ifcW2fc) 

+{Vs - 2Vv){miuj2 + m2UJi)(k^ - mim2 - wifcW2fc)q • k] 

+64(k)[(Vs - 2Vv){uik - 0J2kf<^k^ + Vs(q^ + mim2 + a;ia;2)(k^ - mim2 - a;ifea;2fc)q • k]| . 

(22) 



The normahzation of the component functions is set as 
(fq 16u;iu;2("T'i + ?Ti2)63(q)^4(q) 



(27r)^ (wi + W2)(— q^ + mim2 + ijJ\0J2) 



2M. (23) 



^±/3(q) ^ 7o^±/4(q)i 



2.3 Axial Vector diquark system 

The general form of the wavefunction of 1"^ diquark reads as 

</'(q) = {«±-ei[/i(q)+7o/2(q)+ ^ . ^ 

+M^V5(q) + M;^-^7o/6(q) + #V7(q) + ie^'^'^e^q^jlklM^i)]} 75- (24) 
The constraint conditions are 

a;2 - „ . s , q^ - mim2 - W1W2 , . . , . . ?nia-'2 - m2a;i 

/2(q) = — - — /8(q) + — 777 — - — ^ — /4 q), /e q = ^ttt — - — r/8(q), 

wi + 0J2 M{mi + 1712) M{uii + UJ2) 

M[2Muif5{q) + {1712^1 -miUJ2)fi{q)] ,/ ^ M {miUJ2 + m2Ui) 

= C,^{U,+U2) ' = C,^{U,+U2) 

(25) 

With these conditions the B-S wavefunction of 1+ diquark can be further written as 

q = {«± • ei /i q + 70 — - — /s q + 70 — 777 — - — ^ — /4 q 

0J\ + 002 M{mi + 1712) 
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M[2Ma;i/5(q) + {m2Ui - miu;2)/i(q)] _^ 7o^±/4(q) 



M q2(a;i+a;2) M ^ 

-^x^l ^'^rtT' ^^"" + -"*'^x,7»7.A(q)l}75. (26) 

q (wi + L02) 

The coupled equations are 

n^^^ (^i + '^2)/i (q) - M'^ (mi^a + m2UJi)h (q) 
2(M — — a;2j ; 

OJl +U2 

(q^ + mim2 + '^ia^2)q^/4(q) + M(mi - m2)q^/8(q) , 

Wi + U2 

cfq 



{/i(q)M[(mi - 7712) (Fs + 2Vy){miUJ2k - m2'^ife)(q • k)' 



k2a;ia;2(wifc + ^^2^) 

-(14 - 414)(q^ + mim2 + u;iW2)(wifc + t^2fc)k^q • k] 

+/4(q)k^ [{Vs + 2Vv){miU2 + m2UJi){miuJ2k + m2L0ik)fl ■k-Vs{ui+ U2) {uik + a;2fe) (q • k)^] 

+/5(q)M2[(mi - m2){Vs + 2V^){q^k\uik + u;2k) - 2a;ife(q • k)^) 

+{Vs - 4K)(q^ + mim2 + iOiL02)imiL02k + "i2<^ife)q • k] 

+/8(q)M[Fs (wi +a;2) (m2a;ijk - mia;2fe)q^k^ + (T^ + 2K) (mi(^2 + m2a;i) (wife - a;2fe)k^q • k^ 

(27) 

Mq2(a;i + W2)/i(q) - M^{miU2 + m2Wi)/5(q) 



2(M + a;i + a;2)[ 



uJi + u;2 

(q2 + mim2 + wiu;2)q^/4(q) + M(mi - m2)q^/8(q) 



Ll>1 + CO2 

\s?uJiUJ2{uJik + ^2k) ^~ ^^'^^^^'^^^ ~ "^2)(14 + 2Vv){miU2k - m2<^ifc)(q • k)^ 
-(Vs - 4T4)(q^ + mim2 + uiU2){uik + ^2fe)k^q • k] 
+/4(q)k'[(^. + 2K)(mia;2 + m2U)i){rniU)2k + 'TZ2a;ifc)q • k — 14 (wi + U2) {u)ik + ^2fe) (q • k) ] 

-/5(q)M2[(mi - m2){Vs + 2V,){G^^\^{u:ik + cu2fc) - 2c^u.(q • k)^) 
+ {Vs - 4K,)(q^ + mim2 + ujiuj2){miuj2k + "i'2'^ifc)q ' k] 

+/8(q)-M[ys(a;i +uJ2){m20Jik - miu;2fc)q^k^ + (K + 214) (miW2 + 777-2^1) (wife -u;2fe)k^q- k]}, 

(28) 

M2(q2 + mim2 + a;ia;2)/5(q) - M{ux + a;2)qV8(q) 



4(M — wi - a;2 , 

a;i +U2 

[ V2 d-fii<i)M{Vs + 2Vy){u;i + i02){m2UJik - mia;2fe)[q^k2 - (q • k) 
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+/4(q)k2(mi - m2)Vs{uik + a;2fe)[(q • k)^ - q^k^] 

+2Mk^/8(q)[(ys-2l^^)(q^+mim2+wiW2)(wifc+u;2fc)q-k-V;(mi-m2)(miu;2fc-m2a;ifc)q^] 

(29) 

-M2(q2 + mima + wic^aj/slq) - + W2)qV8(q) 



4(M + a;i +a;2)- 



0^1 + a;2 

k2wiu;2(u)ife + W2fe) 



/ U2 T-^— r{-/i(q)M(y, + 2K)(u;i +a;2)(m2a;ife -mia;2ik)[q'k2 - (q • k)^] 



I (2.)3Mt'i+c.2)2 '^^"^"^^'^"^ + -2)/i(q)/4(q) 



-/4(q)k (mi - m2)Vs{uJik + t^2fc)[(q • k) - q k 
+2M'^f5{ci)[{Vs+2Vy){u;i+u}2){(^2k<i^^^+^^iki(i-^f)-VsimiUJ2+m 

-2Mk^/8(q)[(ys-2F^)(q^+mim2+a;iW2)(wifc+a;2fe)q-k-Fs(mi-m2)(mia;2fe-m2a;ifc)q^]. 

(30) 

The normalization is set as 

(27r)3Af(^i + c^2)2 ■ 
-a;ic^2("iiW2 + m2a;i)M/4(q)/5(q) + u;iL02{miL02 - m2Wi)/i(q)/8(q) 
+M'^{m\ijOi — mim2UJi — mim20J2 + 0Jiu\ — W2q^)/5(q)/8(q)] = 2M. 

3 Numerical results and discussions 

In our numerical computations, the input parameters [21] are iJiu = 0.305 GeV, rrid = 
0.311 GeV, rus = 0.487 GeV, rric = 1.7553 GeV, nib = 5.224 GeV, A = 0.20 GeV^ Kqcd 
0.26 GeV, a = 2.71828 , and a = 0.06 GeV. 

Besides the spectra of different free diquark states, we also evaluate their mean square- 
root radius which is defined as 



V< > = ^<rl> + <rl> + <rl>, (31) 

with 

^ !d^qTr[< y(q)|-^|y(q) >] 

/d3grr[<(^(q)|<^(q)>] " ^ ^ 

Our numerical results are tabulated below. In the following tables, we present the 
results corresponding to different values of (5 and the vacuum expectation value Vq (or the 
zero-point value). 
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VoiGeV) 


0.0 


-0.3 


-0.6 




mass 


msr 


mass 


msr 


mass 


msr 


{ud)o+ 


0.6032 


3.5287 


0.3216 


11.3828 


- 


- 


(n.s)o+ 


0.7839 


3.6089 


0.4962 


11.3616 


0.1982 


11.3802 


(r/.s)o- 


0.7898 


3.6330 


0.5021 


11.3597 


0.2041 


11.3790 


(uc)o+ 


2.057 


5.2610 


1.763 


11.3336 


1.465 


11.3652 


(dc)o+ 


2.063 


5.5278 


1.769 


11.3293 


1.471 


11.3639 


(sc)o+ 


2.230 


2.8423 


1.943 


10.0198 


1.644 


11.2610 


(m5)o+ 


5.527 


5.2038 


5.232 


11.3260 


4.933 


11.3615 


(d6)o+ 


5.532 


6.4324 


5.238 


11.3217 


4.939 


11.3594 


(6s)o+ 


5.698 


2.7435 


5.410 


8.8897 


5.114 


10.8857 


(6c)o+ 


6.918 


1.1725 


6.627 


1.3093 


6.334 


1.4472 


(■uu)i+ 


0.6127 


11.2596 


0.3156 


11.3421 


- 


- 


(wf/)i 


0.G18G 


11.2483 


0.3215 


11.3416 


- 


- 




0.6246 


11.2458 


0.3275 


11.3411 


- 


- 


(iis)i + 


0.7940 


11.2386 


0.4962 


11.2577 


0.1980 


11.0937 


(cis)i+ 


0.8000 


11.0985 


0.5021 


10.5140 


0.2038 


11.2139 




0.9725 


5.5968 


0.6767 


11.2590 


0.3778 


11.2756 


(uc)i+ 


2.060 


5.4201 


1.763 


4.0055 


1.465 


3.8727 


((ic)i+ 


2.066 


5.2305 


1.769 


3.9799 


1.471 


3.8273 


(sc)i+ 


2.232 


2.9041 


1.942 


7.9565 


1.644 


3.5679 


(cc)i+ 


3.460 


1.3865 


3.169 


1.5507 


2.876 


1.7133 




5.528 


4.3950 


5.232 


11.2767 


4.933 


4.8772 




5.534 


4.2068 


5.238 


5.0016 


4.939 


4.7222 


(fe)i+ 


5.700 


2.4584 


5.411 


4.9070 


5.113 


11.1654 




6.912 


1.1162 


6.621 


1.2287 


6.328 


1.3427 


(&fe)i+ 


10.33 


0.7254 


10.04 


0.7500 


9.740 


0.7742 



Table 1: diquark mass(GeV) and mean square-root radius (fm) with (3 = 0.0 
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VoiGeV) 


0.0 


-0.3 


-0.6 




mass 


msr 


mass 


msr 


mass 


msr 


{ud)o+ 


0.7397 


4.0056 


0.6724 


2.0227 


0.5139 


2.3418 


(n.s)o+ 


0.9496 


1.7019 


0.8372 


1.7572 


0.6583 


2.1020 


(d,s)() 


0.9565 


1.6778 


0.8423 


1.7470 


0.6022 


2.0898 


(uc)o+ 


2.274 


1.1711 


2.096 


1.4230 


1.889 


1.7430 


(dc)o+ 


2.280 


1.1652 


2.100 


1.4152 


1.898 


1.7300 


(sc)o+ 


2.445 


1.0478 


2.232 


1.2454 


1.999 


1.4569 


(u6)o+ 


5.759 


1.1265 


5.569 


1.3731 


5.356 


1.6822 


(d6)o+ 


5.764 


1.1206 


5.572 


1.3648 


5.359 


1.6686 


(6s)o+ 


5.920 


0.9992 


5.696 


1.1815 


5.458 


1.3777 


{bc)o+ 


7.088 


0.7089 


6.810 


0.7576 


6.532 


0.8039 


(uu)i+ 


0.8846 


1.5497 


0.6868 


1.6647 


0.4097 


1.4848 




0.8904 


1.5471 


0.6918 


1.6598 


0.4152 


1.1856 


idd)i+ 


0.8962 


1.5376 


0.6968 


1.6538 


0.4206 


1.4857 


{us)i+ 


1.058 


1.4241 


0.8476 


1.5606 


0.5896 


1.5354 


{ds)i+ 


1.064 


1.4190 


0.8521 


1.5538 


0.5935 


1.5294 


{ss)i+ 


1.229 


1.2974 


0.9982 


1.4166 


0.7323 


1.4287 


{uc)i+ 


2.299 


1.2302 


2.097 


1.4450 


1.877 


1.7059 


(dc)i+ 


2.305 


1.2252 


2.101 


1.4375 


1.880 


1.6925 


(sc)i+ 


2.462 


1.0995 


2.228 


1.2476 


1.982 


1.4001 


(cc)i+ 


3.648 


0.8163 


3.371 


0.8633 


3.090 


0.9053 




5.764 


1.1815 


5.567 


1.3683 


5.350 


1.6523 




5.769 


1.1826 


5.570 


1.3887 


5.353 


1.6462 


(fe)i+ 


5.923 


1.0507 


5.694 


1.2232 


5.453 


1.4060 


(bc)i+ 


7.078 


0.7069 


6.798 


0.7468 


6.517 


0.7855 




10.46 


0.5329 


10.17 


0.5455 


9.874 


0.5579 



Table 2: diquark mass(GeV) and mean square root radius (fm) with /3 = 0.25 
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VoiGeV) 


0.0 


-0.3 


-0.6 




mass 


msr 


mass 


msr 


mass 


msr 


{ud)o+ 


0.7904 


1.3035 


0.7911 


2.0982 


0.7105 


1.7727 


(n.s)o+ 


1.028 


2.3651 


0.9731 


1.5194 


0.8582 


1.5598 


(d,s)() 


1.03G 


2.2344 


0.9789 


1.504G 


0.8023 


1.5515 


(uc)o+ 


2.394 


0.9782 


2.251 


1.0897 


2.086 


1.2343 


(dc)o+ 


2.400 


0.9725 


2.255 


1.0838 


2.088 


1.2277 


(sc)o+ 


2.575 


0.8603 


2.393 


0.9666 


2.192 


1.0854 


(u6)o+ 


5.888 


0.9047 


5.728 


1.0192 


5.553 


1.1626 


(d6)o+ 


5.894 


0.8995 


5.733 


1.0130 


5.556 


1.1560 


(6s)o+ 


6.056 


0.7977 


5.859 


0.8979 


5.648 


1.0088 


{bc)o+ 


7.215 


0.5893 


6.946 


0.6250 


6.675 


0.6574 


(uu)i+ 


0.9900 


0.4521 


0.8590 


1.3209 


0.6549 


1.2716 




0.9973 


0.480C 


0.8644 


1.3170 


0.G59C 


1.2697 


idd)i+ 


1.004 


0.5075 


0.8697 


1.3115 


0.6642 


1.2664 


{us)i+ 


1.192 


0.7153 


1.0222 


1.2223 


0.8099 


1.2278 


{ds)i+ 


1.199 


0.7244 


1.027 


1.2185 


0.8139 


1.2239 


{ss)i+ 


1.375 


1.0575 


1.178 


1.1206 


0.9470 


1.1405 


{uc)i+ 


2.428 


0.9021 


2.252 


1.0265 


2.064 


1.1565 


(dc)i+ 


2.434 


0.9030 


2.257 


1.0251 


2.067 


1.1527 


(sc)i+ 


2.603 


0.8764 


2.391 


0.9605 


2.168 


1.0456 


(cc)i+ 


3.789 


0.6819 


3.520 


0.7152 


3.247 


0.7451 




5.884 


0.7958 


5.717 


0.9110 


5.535 


1.0536 




5.890 


0.7996 


5.721 


0.9186 


5.539 


1.0403 


(fe)i+ 


6.056 


0.8078 


5.852 


0.8995 


5.636 


1.0172 


(bc)i+ 


7.206 


0.5962 


6.933 


0.6259 


6.658 


0.6527 




10.56 


0.4669 


10.28 


0.4765 


9.985 


0.4854 



Table 3: diquark mass(GeV) and mean square root radius (fm) with (3 
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VoiGeV) 


0.0 


-0.3 


-0.6 




mass 


msr 


mass 


msr 


mass 


msr 




- 


- 


0.4383 


9.1040 


0.6124 


1.7486 


(ms)o+ 


1.069 


9.1714 


1.059 


1.6103 


0.9856 


1.4060 




1.083 


3.2217 


1.005 


1.5805 


0.9904 


1.3955 


{uc)o+ 


2.485 


0.8955 


2.364 


0.9580 


2.223 


1.0424 




2.492 


0.8894 


2.368 


0.9524 


2.226 


1.0372 


(sc)o+ 


2.677 


0.7734 


2.516 


0.8427 


2.337 


0.9227 


{ub)o+ 


5.987 


0.8029 


5.846 


0.8744 


5.692 


0.9620 


(ri6)o+ 


5.993 


0.7978 


5.850 


0.8692 


5.695 


0.9565 


(6s)o+ 


6.163 


0.7012 


5.984 


0.7688 


5.792 


0.8451 


(6c)o+ 


7.324 


0.5236 


7.063 


0.4804 


6.797 


0.5805 


{uu)i+ 


- 


- 


0.3729 


10.0319 


0.8123 


1.1460 


(ad),. 


- 


- 


0.4387 


9.9951 


0.8172 


1.1435 


idd)i+ 


- 


- 


0.4961 


9.9615 


0.8220 


1.1391 


{us)i+ 


1.069 


10.0515 


1.143 


1.0528 


0.9662 


1.0805 




1.099 


10.0266 


1.149 


1.0528 


0.9707 


1.0776 


{ss)i+ 


1.481 


0.9391 


1.311 


0.9839 


1.108 


1.0013 


{uc)i+ 


2.525 


0.7164 


2.366 


0.8233 


2.195 


0.9183 


{dc)i+ 


2.532 


0.7193 


2.371 


0.8246 


2.199 


0.9182 


(.sc)i+ 


2.714 


0.7469 


2.518 


0.8146 


2.310 


0.8776 


(cc)i+ 


3.909 


0.6086 


3.648 


0.6360 


3.381 


0.6597 




5.975 


0.6235 


5.823 


0.7122 


5.662 


0.7935 




5.980 


0.6426 


5.828 


0.7176 


5.666 


0.7935 




6.158 


0.6628 


5.970 


0.7165 


5.771 


0.8260 


(bc)i+ 


7.315 


0.5350 


7.048 


0.5573 


6.778 


0.5812 




10.66 


0.4299 


10.37 


0.4373 


10.08 


0.4445 



Table 4: diquark mass(GeV) and mean square root radius (fm) with /3 = 0.75 
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VoiGeV) 


0.0 


-0.3 


-0.6 




mass 


msr 


mass 


msr 


mass 


msr 


{ud)o+ 


- 


- 


- 


- 


- 


- 




- 


- 


- 


- 


0.7448 


9.8071 


(d.s)o 1- 


- 


- 


- 


- 


0.7903 


9.8484 


(uc)o+ 


2.561 


0.8500 


2.454 


0.8873 


2.331 


0.9410 


{dc)o+ 


2.568 


0.8435 


2.460 


0.8814 


2.335 


0.9359 


(sc)o+ 


2.762 


0.7208 


2.617 


0.7705 


2.454 


0.8285 


{ub)Q+ 


6.070 


0.7409 


5.942 


0.7918 


5.803 


0.8526 




6.077 


0.7359 


5.946 


0.7869 


5.806 


0.8481 


{bs)o+ 


6.254 


0.6418 


6.087 


0.6919 


5.910 


0.7490 


(6c)o+ 


7.417 


0.4771 


7.153 


0.4883 


6.906 


0.5295 


(uu)i+ 


- 


- 


- 


- 


- 


- 




- 


- 


- 


- 


- 


- 


{dd)i+ 


- 


- 


- 


- 


- 


- 


{us)i+ 


- 


- 


- 


- 


0.7450 


10.5665 


ids)i+ 


- 


- 


- 


- 


0.7905 


10.6108 


{ss)i+ 


0.9872 


10.4853 


1.385 


10.8219 


1.236 


0.9149 


(nc)i+ 


2.606 


0.5984 


2.458 


0.6957 


2.300 


0.7763 


{dc)i+ 


2.613 


0.6016 


2.464 


0.6991 


2.304 


0.7776 


{sc)i+ 


2.807 


0.6531 


2.623 


0.7178 


2.427 


0.7731 


(cc)i+ 


4.015 


0.5604 


3.760 


0.5831 


3.500 


0.6037 


{ub)i+ 


6.048 


0.5339 


5.909 


0.5870 


5.760 


0.6630 




6.056 


0.5290 


5.914 


0.5896 


5.765 


0.6567 


ibs)i+ 


6.244 


0.5782 


6.067 


0.6251 


5.881 


0.6807 




7.414 


0.4878 


7.151 


0.5127 


6.886 


0.5403 


{bb)^+ 


10.75 


0.4052 


10.46 


0.4104 


10.18 


0.4166 



Table 5: diquark mass(GeV) and mean square root radius (fm) with /? = 1.0 



In the approach, we take the parameters which are obtained by fitting data for mesons 
[20, 21], then the only parameters which can be adjusted are /3 and Vq. There are several 
remarks to make. 

(1) The masses of 0"*" and 1"^ diquarks which are composed of two heavy quarks or 
one-hcavy-onc-light quarks tend to be degenerate. This is understandable in the heavy 
quark effective theory (HQET) [22] and the results are consistent with that given in refs 
[4, 23]. 

(2) With the same /3 value, the estimated masses of diquarks increase as Vq is larger. 

(3) The masses of diquarks increase with increment of the /? value. 

It is also observed from the numerical results, that for lighter diquarks the axial diquark 
is 200 MeV heavier than the corresponding scalar diqaurk which is consistent with [24] and 
the lattice calculations [23]. As indicated, our approach might be suspicious for dealing 
with light diquarks, however, consistency of the numerical results with that obtained in 
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other approaches indicates its hmited plausibihty. We will discuss this issue in next section. 



4 Discussion and conclusion 

In this paper wc systematically construct the wavefunctions of the scalar and axial vector 
diquarks and then numerically evaluate their spectra. 

There are several free parameters which cannot be determined so far, because the 
diquarks are not experimentally measurable. If one tries to fix the parameters, he has to 
deal with the diquarks which reside in baryons. However in that case, the diquarks are no 
longer free. As the first step, in our work, we investigate the free diquarks which cannot 
independently exist in real world, and in our later works, we will take into account the 
effects due to existence of the extra quark in baryon. We employ the Bethe-Salpeter theory 
to study their spectra and some characteristics. One of the free parameters is (3 which 
signifies the difference for the non-perturbative QCD confinement effect between diquark 
and meson and from a naive consideration, it should be within a range < /? < 1.0. 
Our numerical results confirm that as /3 > 0.75, there is no solution. Another important 
parameter is Vq which stands as the zero-point energy (in the momentum space, it has a 
form of oc Vb5^(q)). In the case of mesons, we know well that such zero-point energy must 
be introduced to meet the data. In the case of diquark, it plays even more important role 
as it severely determines the characteristics of the two-quark system. 

Moreover, as we apply the B-S approach to evaluate the diquarks of two heavy quarks 
or one-heavy-one-light quarks, everything works well, however, once we extend the same 
approach to the system of two light quarks, the solutions do not seem to be reasonable, 
or there are no solutions at all as the parameters take certain values. It is not surprising, 
because as well known, it is hard to use either the potential model or B-S equations to 
deal with pions which are supposed to contain q and q' with q, q' being u, d quarks. It is 
believed that the quark masses in this case are not simply a constant. The case about 
pions are carefully analyzed by Dai et al. [14] and the diquark system with two light 
quarks are studied by Wang et al. [15]. 

Indeed the diquark is not color-singlet and not a real physical state. It resides, gen- 
erally, in baryons. Therefore, the estimated mass and mean square-root radius of a free 
diquark may be different from its practical value in baryons, because of the QCD interac- 
tion of the extra quark with this subject (diquark), if it indeed exists. That is the goal of 
our next work. 

In the relativistic quark potential model, Ebert et al.[25] obtained the values of the 
masses and < >^/2 3.226 GeV, 0.56 fm and 9.778 GeV, 0.37 fm for l^S-i of the 
axial diquarks cc and bb respectively, with inputs nic = 1.55 GeV and nib = 4.88 GeV. 
Comparing our results with theirs, deducing the rest mass contributions (2mc and 2mb) 
whereas we employ larger input values for nic and rrif, in the numerical computations, our 
spectra are consistent with theirs and the mean square-radius estimated in this work is 
about 1.2 times larger than that given in ref. [25]. Considering theoretical uncertainties, 
our results are qualitatively consistent with the values of ref. [25]. 

For lighter scalar diquarks, many authors estimated their spectra and obtained widely 
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diverging results. For example, the authors of [26, 27, 28] obtained M^^,^)^^ < O.SGeV, 
whereas the authors of [29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39] got M(^^a)^+ = 0.3-0.6Gey, 
but M(„rf)^^ > O.GGeV was suggested in Refs. [40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50]. 
For M(„3)°^, it is estimated as 0.63Gey [29], 0.88Gey [40], 0.948Gey [45] and 0.895GeV 
[51]. For °M^ud),+ , it is estimated as 0.614GeF [31] , 0.806 ~ 0.95GeF [23, 41, 45, 51, 40], 
imGeV ~ 1.27Gel^ [34, 48, 47], > 1.27Gey [52]. For M(„,)^^ , 1.069Ge1^ [45] and 1.5Gey 
[51]. For M(^ss).^+, it is estimated as 1.203GeF[45] and 1.215'GeV [51]. The mass of {cu)o+ 
is given as 1.933GeF [53]. The various numbers indicate large theoretical uncertainties. 
To fix them, one needs to apply the results to baryons which is compose of diquarks and 
are observable physical states. 

As suggested, there could be the 1/2-rule that all the corresponding values and prob- 
ably Vq employed for the diquarks, i.e. f3 and Vq should be 1/2 of the values for mesons. 
This rule follows the assumptions that the governing interaction is QCD, thus both short- 
distance (the Coulomb) and the long-distance (the confinement) are proportional to the 
expectation value of the Casimir operator < A°A" > where A" is the SU(3) Gell-Mann 
matrix. If so, the set of results with /3 = 0.5 should be more reasonable. However, the 
situation may be more complicated. Thus in this work, we keep it as an adjustable free 
parameter. As we indicated above, the diquark is not a physical state, so that to fix the 
parameter we need to put them into the baryons, or may be in the future, as hoped, in 
the high-energy heavy ion collisions, free diquark might be directly produced, and then 
we can have more information about its structure. 

Indeed, until we know how the diquark interacts as a whole object with gluons, we 
cannot more reliably estimate the baryon case and evaluate the production rate and decay 
width of the baryon in the quark-diquark picture. Therefore, in our coming work, we are 
going to derive all the form factors at the effective interaction vertices as the diquarks are 
treated as a whole object. 

Acknowledgement: This work is partly supported by the National Natural Science Foun- 
dation of China (NNSFC) and The Special Fund for the Ph.D programs of Chinese Univer- 
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